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Abstract 

We discuss the covariant formulation of the dynamics of particles with abelian and non-abelian gauge charges 
in external fields. Using this formulation we develop an algorithm for the construction of constants of motion, 
which makes use of a generalization of the concept of Killing vectors and tensors in differential geometry. We 
apply the formalism to the motion of classical charges in abelian and non-abelian monopole fields. 
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1 Introduction 



A standard problem in the solution of classical and quantum mechanical systems is to identify 
the constants of motion associated with the system. There are many ways to obtain constants 
of motion, but two methods are fairly common in different branches of physics. In the 
canonical phase-space formulation of classical conservative systems we have to identify all 
dynamical quantities Q which commute with the hamiltonian in the sense of Poisson brackets^: 

{Q,H} = 0. (1) 

The draw-back of this prescription is, that for systems with gauge interactions this formulation 
is usually not manifestly gauge covariant. 

On the other hand for systems with a non-flat configuration space, such as particles 
moving on a curved manifold (or space-time, in general relativity) the appropriate algorithm 
is to search for Killing vectors and their higher-rank generalizations. In Riemannian geometry 
these are covariant objects, but the procedure is only applicable for geodesic motion in the 
absence of non-geometrical external fields of force. 

As observed in pQ, for constants of motion to exist in the case of non-geodesic motion, 
e.g., for particles in external fields, the symmetries of the metric and those of the external 
fields have to match. In fact Killing vectors appear explictly in the expressions for constants 
of motion linear in the momentum. In [21 [3] a complete set of consistency conditions for the 
existence of constants of motion were derived for particles in arbitrary background geometries, 
using a covariant hamiltonian phase-space approach including the contributions of spin. This 
procedure also applied to constants of motion which are higher-order polynomials in the 
momentum, as well as constants of motion which are Grassmann-odd expressions in the 
spin degrees of freedom, which generate standard or non-standard supersymmetries on the 
worldline of the system. 

In this paper we show how to extend this covariant phase-space approach to include 
the presence of external gauge fields. As in [1\ our non-abelian dynamics is based on the 
equations of motion postulated by Wong [5J. These equations were studied in a geometric 
setting, using the method of co-adjoint orbits, for a similar purpose in [5], whilst a lagrangean 
realization in terms of Grassmann variables was constructed in [6]. Having a completely 
covariant phase-space formulation we derive a set of generalized Killing equations, the solution 
of which produces all constants of motion in a manifestly covariant way. To avoid unnecessary 
complications, we formulate all our dynamical models in euclidean or riemannian space, but 
the generalization to minkowskian or lorentzian manifolds is straightforward. 



2 Dynamics of point charges 

The classical dynamics of a point charge in a magnetic field is described by the Lorentz force 
law 

mv = gvxB (2) 
In the standard canonical formulation this equation is derived from a hamiltonian 



1 In quantum mechanics, the Poisson bracket is replaced by the commutator; for brevity and to avoid 
irrelevant operator ordering complications, we stay with classical mechanics in this paper. 
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via Hamilton's equations: 

r= — = -(p- g A), P = — s- = - (VA) • (p - gA) (4) 
ap m or m 

Therefore 

p = mr + gA = mv + gA, p = gVA-v, (5) 

and after substitution and the definition B = V x A eq. ([2]) follows. In terms of the field- 
strength tensor F 

Fij = e ijk B k = ViAj - VjAi, (6) 
the equation for the Lorentz force takes the form 

mvi = qFijVj. (7) 

In this form the equation can be extended easily to relativistic particles in Minkowski space. 

The above construction uses cartesian co-ordinates r*j and their canonical momenta pi, 
such that the equations of motion an be written in terms of Poisson brackets 

' dri dpi dpi dri ' 

i 

for phase-space functions f(j"i,Pi,t) and g(ri,pi,t). In term of these brackets 

The phase-space co-ordinates (r{,Pi) are canonical as the only non-trivial fundamental bracket 
is 

{ri,Pj} = Sij, (10) 

all others vanishing: 

{n,rj} = {pi,Pj} = 0. (11) 
A disadvantage of this formulation is, that the canonical momenta are gauge dependent: 

A' = A + VA => p' = p + gVA, (12) 

although this does not affect the fundamental brackets (I10p . As a result, the hamiltoninan 
equations of motion are not manifestly gauge covariant. 

However, an alternative exists in which the dynamical variables of the particle are all 
gauge invariant, and which has the added advantage that the hamiltonian takes a very simple 
form. Introduce the gauge-invariant momenta 

II = p — qA = mv. (13) 

Then the hamiltonian takes the simple quadratic form 

h = — n 2 . (14) 

2m v ; 

This has the form of a free-particle hamiltonian, but the dynamics is now manifest in the 
modified brackets: 

, df dg df dg df dg 

{f > 9} = ^dfu ~ dfuW t + ^mdny (15) 
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In particular the fundamental brackets are 

{r i ,U j } = S ij , { n ,r i } = 0, {U i ,U j } = -qF ij . (16) 

This shows, that the momenta II are not canonical, but act like covariant derivatives, rather 
than ordinary partial derivatives; indeed, the last bracket is the Poisson-bracket version of the 
Ricci identity. As a result, we can derive the homogeneous Maxwell equations (the Bianchi 
identities) from the Jacobi identity: 

{U i ,{U j ,U k }} + {U j ,{U k ,U i }} + {U k ,{U i ,U j }} = (17) 

which implies 

ViFjk + VjFki + VkFij = & V-B = 0. (18) 

It remains to establish that the brackets and the hamiltonian reproduce the correct equations 
of motion; this follows by direct computation: 

h = {r u H} = 5i, fli = {ILi, H} = ± FijILj. (19) 
m m 

The covariant phase-space formulation has been used by various authors in different contexts, 
see e.g. [HESE]. 



3 Symmetries and constants of motion 

The gauge-covariant formulation of hamiltonian mechanics of charged particles is mathemat- 
ically elegant; it is also most suited to study the existence of symmetries and constants of 
motion. Indeed, in the hamiltonian framework a constant of motion Q(r, II) is identified by 
the property that its bracket with the hamiltonian vanishes: 

{Q,H} = => Ui(ViQ - qFi, = n • (VQ + qB x = 0. (20) 

A systematic procedure is to expand any constant of motion as a power series in II: 

Q(r,n) = Civ) + C i (r)II i + ^0^)11^ + ... . (21) 

Substitution gives a series of constraints 
VjC = qFijCj, 

ViCj + VjCi = qFikCkj + qFjkCki, 

(22) 

VjCjfc + VjCki + ^kCij = qFuCijk + qFjiQki + qFkiQij, 



This series can be truncated whenever there is a Killing vector or tensor of flat space: 

V (ii q,.., n )=0. (23) 
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Then we can take C^.; = for all p > n, and the constant of motion takes the polynomial 
form 

P-l - 

Q(r,U) = ^ - C n ... tk (r)U n ...U tk . (24) 

k=0 

Note that it is always possible to add an arbitary constant to the zeroth order coefficient 
C (r) . Therefore it is obvious that for particles in an electromagnetic background there are no 
non-trivial constants of motion corresponding to only a C(r) with Cj(r), Cj,-(r) and all higher 
co-efficients vanishing. The first non-trivial case is therefore the set truncated at p = 2: 

Q(r,n) = C(r) + C(r) • II, (25) 

with C a Killing vector of flat space: 

ViCj + VjCi = 0. (26) 

Such Killing vectors generate translations and rotations, and take the form 

C = m + nxr, (27) 

where m and n are arbitrary constant vectors. The lowest-order constraint equation now 
becomes 

ViC = qFijCj = qVi {AjCj) - qAjViCj - qC 3 V 3 Ai (28) 
Now use the Killing equation for Cj to rewrite this equation as 

Vi(C-gC- A) =qA- VC { - qC ■ V A t . (29) 

A very simple example a constant magnetic field: 

A = |C => B = (3n. (30) 



For such a field 



C = gC-A = yC 2 (mod constant). (31) 



As a result the full constant of motion takes the form 

Q = C-(gA + n) = C-p 

= m'p + ^B'(rxp). 

P 



(32) 



As m is arbitrary, all components of the momentum and the component of the angular 
momentum in the direction of B are conserved. This reflects the invariance under translations 
and transverse rotations in a constant magnetic field. 

Starting from the general Killing vector (|27p the equation for the Killing scalar becomes 

VC = qui x B + g(rn • B - nr • B) . (33) 

As a non-trivial example, consider axially symmetric fields B with the axis defined by a unit 
vector n: 

B _ ^(p) r x n = ^(p) x n p = r -(r n)n, (34) 
P P 
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where p = |p|. Then equation ([33]) becomes 

qB 

VC = gm x B = (p n • m — n p • m) . (35) 

P 

The corresponding vector potential can be taken as 

A = g(p)n, (36) 
provided we identify B{p) = g'(p). Then eq. (|35[) becomes 

VC = — — (p n • m — n p • m) = q (m • n) Vg — q (m • Vg) n. (37) 
p 

This equation allows solutions for m = An, with A an arbitrary scalar parameter. As n p = 
and n = 1, one finds 

C = Xqg(p). (38) 

The full constant of motion then reads 

Q = Xqg(p) + (An + n x r) -II 

(39) 

= An • p + n • (r x p) . 

As A is arbitrary, it follows that the components of the canonical momentum and the angular 
momentum in the direction n are independently conserved^! 

One can also search for constants of motion which are higher-order polynomials in the 
momentum, associated with flat-space Killing tensors. The simplest one is 

Cj <%, (40) 

which has the special property that 

FijCjk + FkjCji = 0. (41) 

Therefore the associated Killing vector and scalar can be taken to vanish, and the correspond- 
ing constant of motion is the hamiltonian: 

C = i dijUiUj = mH, (42) 
More complicated Killing tensors are of the form 

Cij = 2 4,-n • r - (rtiTj + njn), (43) 
for an arbitrary fixed unit vector n, and 

= %r 2 -nrj, (44) 

which is the radial counterpart of (03]). Any constants of motion associated with these Killing 
tensors are extensions of the n-component of the Runge-Lenz vector 

nK = n(nxL)=nrII 2 -nIIrn, (45) 

or the total angular momentum 

L 2 = (r x n) 2 = r 2 n 2 - (r -n) 2 = r • K. (46) 

Such constants of motion are associated with special field configurations, in particular spher- 
ically symmetric ones. We discuss such fields in the next section. 

2 Actually, the components of the canonical and covariant angular momentum in the direction n are the 
same. 
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4 Magnetic monopoles 

A spherically symmetric magnetic solution of the Maxwell equations is the Dirac monopole: 

B = (47) 

For such a field equation (|33p takes the form 

VC = ±| (m x r + (r • n)r - r 2 n) . (48) 

Now the first term is a curl, not a gradient; as a result we have to take m = 0, and 

C = -qg—. (49) 
r 

The result for the constants of motion based on the Killing vector (|27|) then is 

g = n-(- 9 ^+rxn), (50) 

where n is an arbitrary vector; therefore all components of the gauge-covariant improved 
angular momentum 

r 

J = L — grx A - gq-, L = r x p, (51) 
r 

are conserved. We observe, that these quantities generate rotations and satisfy the standard 
so(3) Lie algebra 

{Ji,Jj} = EijkJk- (52) 
The Casimir invariant of so(3) is the total angular momentum squared: 

J 2 = L 2 + 5 V, (53) 

which is a constant of motion with 

C iS = 8iji*-Tirj, d = 0, C = g 2 q 2 . (54) 

It follows from (|53p . that the values of the total angular momentum in the classical theory 
satisfy the bound J 2 > C = g 2 q 2 ■ A recent analysis relating many different formulations of 
this dynamicalsystem is found in [7j. 

Remarkably, the Runge-Lenz vector can be extended to a constant of motion in another 
type of central magnetic field: 

Indeed, in such a field there is a constant of motion 

Q = n.(K + ^L-2gV r -), (56) 



for any unit vector n, with 



Cij = 2 Sij n • r - (mrj + nju), 

7f ' ( 57 ) 

C = -2g 2 q 2 ^. 
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In contrast, for such a central field there is no conserved angular momentum vector J, al- 
though the total angular momentum J 2 is conserved. Of course, the total energy of the 
magnetic field diverges both at r = and at r — > oo, and the field does not satisfy the 
free Maxwell equations. Therefore it requires non-trivial magnetic sources and boundary con- 
ditions. However, even if the field would exist only in a restricted region of space, the constant 
of motion exists provided the orbit of the point charge is also restricted to this region. 

5 Non-abelian point charges 

The gauge-covariant dynamics of point charges can be extended to non-abelian point charges. 
The starting point is defined by the Wong equations [3], which can be written in the form 

II = mv, 11 = gt a v x B a , 

(58) 

ta = gfabcV ■ A b t c , 

where the t a are the non-abelian gauge variables, and where the non-abelian field strength is 
defined as 

Fija = Eijk Bj~a = ^iAj a — VjAj a — gf abc A^ b Aj c . (59) 

A lagrangean representation of the gauge variables in terms of Grassmann-odd degrees of free- 
dom was shown to exist in ref. [6]. The hamiltonian formulation uses canonical co-ordinates 
and momenta, with a hamiltonian 

H = ^- (p - gA a t a ) 2 , (60) 
2m 

supplemented by the fundamental brackets 

i r i,Pj} = {t a ,t b } = -fabctc, (61) 

with all other brackets vanishing. It is easily verified, that these brackets reproduce the 
equations of motion (|58p . A more transparent and gauge-covariant formulation is obtained 
by introduction of the covariant momentum 

II = mv = p - gA a t a . (62) 

with the quasi-free hamiltonian 

H = — n 2 . (63) 

2m v ; 

The equations of motion ([58]) are now reobtained from the covariant brackets 

f = {f,H}, (64) 

where the brackets are defined explicitly by 

, dh df , ^ df dh df dh . N 

{/ ' h} = Vif ml ~ m 1 + 9 ija a mWr fabctc dTadf b (65) 

Here the covariant phase-space derivative of a function /(r, II, t a ) appearing on the right-hand 
side is defined as 

df 

Vif = Vif - gf abc t a A ib —. (66) 

ot c 
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As in the abelian case we can now look for constants of motion by applying the hamiltonian 
formalism: 

dQ 

or in vector notation 



{Q, H} = Ui[ViQ- gF l3 a t a ) = 0, (67) 



n • (VQ - gt a f abc A b ^- + gt a B a X ^) = 0. (68) 
Using a covariant momentum expansion 

Q = C + dlU + \ Cijlkllj + .., (69) 

we obtain a set of constraints to be satisfied 
T^iC = gt a Fij a Cj , 

T>iCj + VjCi = fft a {FikaCkj + gFjkaCki) , 

DiCjk + VjCki + V k Cij = gt a (FnaCijk + Fj ia Ciki + F k i a Cuj) , 



(70) 



Clearly the condition for the series (|69p to stop after a finite number of terms is the existence 
of a tensor satisfying the condition 

%C i2 ... in) =0, (71) 

a gauge-covariant generalization of the Killing equation. All of this is a direct non-abelian 
generalization of the case of Maxwell-Lorentz theory. 

6 2-D Yang-Mills theory 

A simple example to illustrate the procedure of constructing constants of motion from gener- 
alized Killing-vectors and tensors is offered by 2-D SU (2) Yang-Mills theory. In this theory 
the magnetic field strength is represented by a triplet of scalar fields B a : 

F?j = (72) 

which satisfies the euclidean Yang-Mills equations 

D t B a = 0. (73) 

Thus the magnetic field is covariantly constant, and its modulus B 2 = B a B a is a gauge- 
invariant real constant: 

V l B 2 = 0. (74) 

As the effect of a local gauge transformation is a local rotation of B a in the internal space, it 
is possible to gauge transform the magnetic field locally into a constant, e.g. 

B a = (0,0,B). (75) 
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Such a constant magnetic field can be constructed from the linear gauge potential 

A? = ~B a e l3 r r (76) 

In a constant magnetic field there is translation and rotation invariance, hence we can look 
for Killing vectors of the type (|27|) : 

d = mi — Xeijrj, (77) 

with rrii and A arbitrary constants. It then only remains to solve for the generalized Killing 
scalar C = C a t a : 

{VtCf = ViC a - ge abc A\C c 

(78) 

= gF^Cj = gB a (eijmj + An) . 
The straightforward solution to this equation with Af given by (|76p is 



C a = gB a {sijnmj + ^ r 2 J . (79) 

With rrii and A we thus associate two constants of motion; a gauge-improved momentum: 

Qi = n { - gB a t a e ijrj = Pi + gA?t a , (80) 
and the canonical angular momentum 

J = EijnUj + I B a t a r 2 = Eijnpj . (81) 
As a next step, one can look for symmetric Killing tensors C^-; however, the two candidates 

Cij = Cij = 5ijY — TiTj, (82) 

lead us back to the hamiltonian (|63p and the square of the angular momentum J 2 , respectively. 
A non-trivial constant of motion of this type is a Runge-Lenz-like vector 

Ki = r.tf-n.n.ri gBH a Q e^r 2 + r i£jk r k U k ) + X - (gB a t a ) 2 r< r 2 

(83) 

= n p 2 - Pi p • r + i gB a t a Jri-^ (gB a t a f r t r 2 . 

These constants of motion, described simultaneously by the arbitrary linear combination n-K, 
are constructed in terms of the generalized Killing tensors 

Cij = 25ij r • n — ritij — VjUi, 

d = -gB a t a Sij (r jr ■ n + i n jr 2 ^ , ( 84 ) 
C=\ (gB a t a ) 2 r 2 r-n. 
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This solution of the 2-D Yang-Mills equations can be embedded straightforwardly in 3-D 
Yang-Mills theory by taking connections 

A« = - l -B a y, Al= l -B a x, (85) 

whilst all other components of vanish. Actually this amounts simply to embedding an 
abelian solution of the type (|30p into a non-abelian model. In essence, the results of the 
abelian case are reproduced. 



7 The non-abelian Wu-Yang monopole 

A genuinely non-abelian static solution of the pure SU{2) Yang-Mills equations in 3-D space 
is the non-abelian Wu-Yang monopole [Hj; the monopole field is given by 

A a = ISiakTk^ (86) 

g r l 

with the corresponding magnetic field strength 

*g = lS2jl£, D jF ^ = 0. (87) 

In such a field a particle has a conserved charge invariant under combined spatial and isospin 
rotations 

r a t 

Q = — => V l Q = 0. 

r 

In addition, there is a conserved angular momentum 

n • J = n • ( r x II - Q- ) 44> J = rxp-t = L-t. 



r , 

which is constructed from the Killing vector and associated scalar 

a = (nxr)„ C = -C— = -n.r-^. (90) 

It is easily established that the components of the angular momentum J generate the so(3) 
Lie algebra (|52p . The contribution of isospin to the orbital angular momentum mixes gauge 
and spin degrees of freedom, a result well-known in the literature [10] . For point-particles 
carrying even-dimensional isospin representations (e.g., isodoublets) this turns the bound 
states into fermions [11] . 

Also in this case there exist constants of motion quadratic in the momenta: the hamilto- 
nian H and the square of the total angular momentum: 

j 2 = r 2 n 2 - (r • n) 2 + {~~\ 2 ■ (91) 



This constant of motion is constructed from the Killing tensor 

with d = and C = Q 2 . In contrast, a constant of motion polynomial in the momenta 
which generalizes the Runge-Lenz vector does not exist for a point-particle in the Wu-Yang 
monopole background. 

As is well-known, there also exist abelian monopole solutions in spontaneously broken 
non-abelian gauge theories |14} I15j . The motion of non-abelian point particles in such a 
background has been studied in [T7]. 
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8 Charged particles in curved space 

The concept of Killing vector has its origin in differential geometry, where it arises as generator 
of an isometry. In the previous sections we have applied the concept in flat space, the 
isometries of which are translations and rotations. We have shown how the concept can 
be generalized in the presence of background gauge fields, abelian as well as non-abelian. 
Symmetries and constants of motion arise in particular when the isometries are matched by 
symmetries of the background fields. We have also discussed Killing tensors of higher rank, 
associated with constants of motion depending on higher powers of the momenta. 

The generalizations can easily be extended to non-flat spaces. The hamiltonian of a 
charged particle moving in a space with metric gij (x) is 

H = ^9 ij (x)U i U j . (92) 

For a particle without spin the covariant brackets are the same as in flat space: 

{x i ,U j }=S}, {U i ,U j }=qF ij . (93) 
In particular, the equations of motion become 

x i = {x\ H) = — g ij Uj, Ui = mgij x j , 

(94) 

ri = {IT , H] = 1 g kl (iyiljllj + qF ik U^ I ^±=Tl i -x k vJlLj = qF tj j? . 

By the first of these equations, the second one reduces to the Lorentz-Wong equations in 
curved space. As before, constants of motion are obtained by solving the equation 

{Q,H} = 0, (95) 

with Q a polynomial in the momenta 

Q(x,U) = C(x) + C\x)U t + X - C^'(x)ILIL- + ... (96) 

Then the coefficients are solutions of the hierarchy of differential equations 
D i C = V i C = qF ij Ci, 
D t Cj + Djd = qF lk C k 3 + qF ]k C% 

(97) 

DiCjk + DjC k i + D k Cij = qFuC 1 - k + qFjiC 1 ki + qF ki C l {j , 



As usual, indices are raised and lowered with the metric, and the covariant derivative Di is 
constructed with the Levi-Civita connection 

DiCj = ViCj - r./Cfc, (98) 
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in the case of abelian background gauge fields. In the case of non-abelian background gauge 
fields we have to make the replacements 

d 

Di -> T>i = Di - gf abc t a A ib —, qFij -> gt a Fi ja . (99) 

As an example we consider the motion of a charged particle on the unit sphere S 2 supplied 
with a constant magnetic field. The metric on the sphere is defined by the line element 

ds 2 = g {j dx l dx j = d9 2 + sin 2 9 dip 2 . (100) 

The sphere S 2 admits a triplet of Killing vectors C % = (C 6 \C V ) satisfying the homogeneous 
form of the second Killing equation (|97p . with and all higher Killing tensors vanishing: 

C/j\ = (— simp, — cot 9 cos ip), 

CLx = (cos ip, — cot 9 simp), (101) 

q 3) = (°> i)- 

These Killing vectors generate three independent rotations on S 2 . The magnetic field with 
constant flux B is described by the field strength 

Fij dx l A dx j = B sin 9 d9 A cfy>. (102) 

Applying the vectors (|10ip in the right-hand side of the equation for the Killing scalars, we 
find that they take the form 

Cm = — qB sin 9 cos p, Ci 2 \ = — qBsm9simp, Cy 3 \ = — qB cos 9. (103) 

Therefore we find as constants of motion the components of the gauge-improved angular 
momentum Ji a )'- 

Jm = — sin ip Hg — cot 9 cos 93 11^ — qB sin 9 cos ip, 

J(2) = cos ip He — cot 9 sin cp 11^ — qB sin 9 sin <p>, (104) 

J( 3 ) = - qB cos 9. 
As in eq. (|52p . these constants of motion satisfy the so(3) Poisson bracket algebra 

{J(a), J(b)} = £ abc J(c)- (105) 

Indeed, the present model is equivalent to a dimensional reduction of the monopole field from 
3-D flat space to the 2-D unit sphere [7]. As a result, we also expect the existence of Killing 
tensors. First observe, that 

J(a) = 2mH + <^ 2 , (106) 

a 

which is a quadratic expression in the momenta with Cij = gij. In addition, there two other 
independent symmetric Killing tensors: 

c ij = f cos ip \ clj = f s\mp \ 

( x ) y cosy? —2 sin ip cot 9 J ' ( 2 ) \ sin ip 2 cos ip cot 9 J 
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Inserting these expressions on the right-hand side of the second equation (|97p . and solving 
this equation, we find the associated generalized Killing vectors 



qB 

C/i \ = (cos 9 sin 9 cos p, — (cos 2 9 — sin 2 9) sin p) 

' sin 9 

qB 

CL^ = (cos 9 sin 9 sin p, (cos 2 9 — sin 2 9) cos p) . 

1 ; sin9 11 



(110) 



(108) 



The associated Killing scalars, the solution of the first equation (|97p . read 

C(X) = q 2 B 2 sin 9 cos 9 sin p, CV 2 ) = — q 2 B 2 sin cos cos 99. (109) 
Combining these results we find the constants of motion 

2 n cos sin 99 
K(d = cos pUgllu, — cot 9 sin pTl. n — qB cos 9 cos pile + qB : — H v 

+ g 2 i? 2 sin cos sin </>, 

2 /1 cos 2^ cos p 

K( 2 \ = sin 03 rhjllo + cot 9 cos pH w — qB cos 9 sin pTie — qB H v 

w * sm9 

— q 2 B 2 sin # cos 9 cos p. 

Observe, that 

K {1) = d v K {2) , K {2) = -8 V K (1) . (Ill) 
These relations follow, because d^Di = Did v . 

9 Super symmetry 

Spinning particles whose internal angular momentum is described by Grassmann co-ordinates 
if) 1 can have Grassmann-odd constants of motion, generating transformations in anti-commuting 
co-ordinates. If their bracket closes on the hamiltonian, they generate standard supersym- 
metries. In the case of charged particles in an external gauge field, the standard supercharge 
takes the form 

n = n t ip\ (112) 

whilst the non-zero covariant brackets are 

{x i ,U J }=5}, {U i ,Il j } = qF ij , {^\^} = -i8 i K (113) 

in the abelian case, with appropriate modifications in the non-abelian generalization. It 
follows, that the internal spin rotations are generated by the bilinears 

Si = £ijk tp j ip h , {si, sj} = e ijk s k . (114) 



The hamiltonian in flat space reads 



H = — IT 2 - — B s, (115) 
2m m 
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and satisfies the supersymmetric bracket relation 

{n,n} = -2miH. (116) 

In such a theory any dynamical quantity of which the bracket with the supercharge vanishes, 
is automatically a constant of motion: 

{Q,n} = {Q,H} = 0, (117) 

owing to the Jacobi identity for the brackets. The reverse does not hold in general. Hence 
the class of super invariants is a subclass of the constants of motion. For these superinvariants 
one can derive another more restrictive hierarchy of conditions which are sufficient, though in 
general not necessary, to obtain solutions of equations (f22|) or their appropriate generalizations 
(fTOj) or (|U7|) . These equations were derived in [3], hence it is not necessary to elaborate on 
them in detail. The generating equation is 

-^(v,C?- 9 F,||) + n,§|=0, (118) 

obtained by writing out the bracket {£l,Q}. The hierarchy of square roots of the extended 
Killing equations is obtained by expanding Q in a series in the momenta IT, and each coef- 
ficient Ci 1 ...i n (x,vp) in a (finite) polynomial in the Grassmann variables ip l . In some cases of 
physical interest the superinvariants do not only include known constants of motion, such as 
the angular momentum in the case of a monopole field, but also new Grassmann-odd invari- 
ants. The coefficients in the expressions for such conserved odd charges are generalizations 
of the so-called Killing- Yano tensors, rather than the Killing tensors proper. In the case 
of the magnetic monopole such an anti-commuting constant of motion is the non-standard 
supercharge p2] 

Q = e ijk (xWi/j* - l - ^W") , (H9) 

the bracket of which with itself closes on the square of the angular momentum, rather than 
on the hamiltonian: 

{Q,Q}=-i(j 2 -2gq r -^), (120) 
with J = r x II + s. A generalization for the non-abelian monopole was constructed in |13j . 

10 Summary and discussion 

In this paper we have developed an algorithm to construct all constants of motion for con- 
servative dynamical systems. The method, based on extensions of the Killing equations in 
differential geometry, works in the presence of gauge interactions and in non-flat geometries 
as well. It brings out in particular the importance of tuning the symmetries of the external 
fields with those of the geometry of the configuration space. The method has been illustrated 
with several examples, in particular monopole-type solutions in abelian and non-abelian gauge 
theories. We have restricted ourselves to classical dynamical systems, but the use of a bracket 
formulation on phase space allows easy translation -modulo operator ordering- to the case 
of quantum systems. Also, we have not included spin degrees of freedom. Extending the par- 
ticle models with Grassmann variables to describe fermions opens the possibility to include 
super symmetries in this framework. 
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Apart from the generic importance of constructing constants of motion, the classifica- 
tion of all dynamical variables commuting with the hamiltonian is a starting point for the 
procedure of hamiltonian reduction. This procedure provides an elegant way of constructing 
non-trivial integrable models; for a recent review see [18]. This technique was applied to derive 
N = 4 supersymmetric mechanics in a monopole background in [TP]. The connection with 
Killing vectors and tensors and their generalizations discussed here provides a geometrical 
configuration-space description of this reduction process. 
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